Configuration of separability and tests for multipartite entanglement in Bell-type 

experiments 
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We derive tight quadratic inequalities for all kinds of hybrid separable-inseparable n-particle den- 
sity operators on an arbitrary dimensional space. This methodology enables us to truly derive a 
tight quadratic inequality as tests for full n-partite entanglement in various Bell-type correlation 
experiments on the systems that may not be identified as a collection of qubits, e.g., those involving 
photons measured by incomplete detectors. It is also proved that when the two measured observables 
are assumed to precisely anti-commute, a stronger quadratic inequality can be used as a witness of 
full n-partite entanglement. 



PACS numbers: 03.65.Ud 

Sincc 1980's, it has been a problem how to confirm mul- 
tipartite entanglement experimcntally. Recently, we have 
been given precious experimental data by efforts of ex- 
perimentalists |ï],|| . Proper analysis of these experimen- 
tal data then becomes necessary, and as a result of such 
analysis the experimental data obtained by Pan and 
co-workers Q confirms the existence of genuinely three- 
particle entanglement under the assumption that proper 
observables are measured in the experiment. However, it 
was discussed [Q that for other experimental data there 
is a loophole problem in confirming three-particle entan- 
glement, and the loophole problem remains unresolved. 
This means that there have not been enough discussions 
about what kind of data are needed for confirming mul- 
tipartite entanglement. 

There have been many researches on the problem that 
provide incmialitics for functions of experimental correla- 
tions Pflp|Jqj7|,p|,p|,^0[. Among them, assuming fc-partite 
split of the system without assuming a specific parti- 
tion, Werner and Wolf derived an upper bound 2^ n ^ k ^ 2 
for expectation vàlues of n-particle Bell-Mermin opera- 
tors B, B' under the assumption that suitable par - 
tial transposes of the density operator are positive ppj . 
The inequality derived by Werner et al. is useful because 
it telis us a number k such that the given state is at most 
fc-separable |11 Ï3]. 

Recently, Uffink has introduced a non-linear inequality 
which is aimed at giving stronger tests for full n-partite 
entanglement than previous formulas. For qubit systems, 
Uffink has derived || a tight quadratic inequality for the 
states where one qubit is not entangled with any other 
qubit; namely, the states written as a convex sum over 
the states of form pi <g> p2,..., n - 

In most of real experiments, we have to deal with 
higher dimensional systems rather than qubit systems. 
For example, when polarizations of photons from a non- 
ideal source are measured by imperfect detectors, it is 
difficult to claim strictly that the observed correlations 
are obtained by measuring subsystems with only two- 
dimensional Hilbcrt spaces, due to the ambiguity in the 
number of photons. The arguments about higher dimen- 



sional systems will thus be necessary in order to establish 
tests applicable to real experiments without making aux- 
iliary assumptions as to the dimension of the measured 
space or as to measured observables. 

In deriving a witness of full n-partite entangle- 
ment, it should be ensured that the witness rules 
out all hybrid separable-inseparable states except gen- 
uine fully n-partite entangled states. The hybrid 
separable-inseparable states are depicted as follows. Con- 
sider a partition of n-particle system {l,2,...,n} into 
k nonempty and disjoint subsets a%, . . . jOt/., where 
Sí=i \ a i\ ~ n J to which we refer as a A;-partite split of the 
system . Let us now consider the density operators p 
on TL = where Hj represents the Hilbert space 

with respect to particle j. Then all hybrid separable- 
inseparable states with respect to partition a\,...,ak 
can be written as 



(d 



where p" 1 , VZ are the density operators on the partial 
Hilbert space ® JgQ . Hj. States (0) are called fc-separable 
with respect to a partition a±, . . . , otk- 

In this paper, we derive the optimal upper bound of 
(B) 2 + (B'} 2 for any partition of the systems ai, . . . , a* 
of an arbitrary dimensional space. It turns out that the 
optimal upper bound depends only on two parameters 
k and rn, and not on the detailed configuration of the 
partition, where m is the number of partides which are 
not entangled with any other partides. The maximum is 
given by 2™+ m ~ 2fe+1 . Using this maximum, we genuinely 
prové that the optimal upper bound that is utilizable to 
confirm full n-partite entanglement (n > 3) of an arbi- 
trary dimensional system is 2 11-2 . Later, we show that 
if an auxiliary assumption as to measured observables is 
allowed, a stronger quadratic inequality can be used as a 
witness of full n-partite entanglement. 

In what follows, we derive tight quadratic inequalities 
for hybrid separable-inseparable states with respect to 
partition ex\, . . . ,otk of an arbitrary dimensional space. 
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It is assumed that a measurement with two outcomes, 
±1, is performed on each particle. Such a measurement 
is generally described by a positive-operator-valued mea- 
sure (POVM), {F+, F+ + F_ = 1, F+, F_ > 0, and 
the corresponding observable is given by a Hcrmitian op- 
erator A = F + — F-, which has a spectrum in [—1,1]. We 
assume that for each particle j, either of two such observ- 
ables Aj or A'j is chosen, where — 1 < Aj,A'j < l,Vj. 

The Bcll-Mcrmin operators take a simple form when 
we view on a complex plane using a function f(x,y) = 
-^e~ 17r / 4 (.T + iy),x,y £ R. Note that this function is 
invertible, as x = 3?/ - y = 3?/ + 3/. The Bcll- 
Mcrmin operators £?n„ and B' N arc dcfincd by |ïo| , [Ï2"f 

where N„ = {1, 2, . . . , n}. We also define B a for any 
subset a C N n by 

f(B a ,B' a ) = ® jea f(A j ,A' j ). (3) 

It is easy to see, when a,/J(c N n ) are disjoint, that 

f(B aU0 ,B' aup ) = f(B a ,B' a ) ® f(B p ,B' p ), (4) 

which leads to 



B, 



aU/3 



1/2(B Q B^ + B'jBp) + l/2{B a Bp - B' a B' p ), 



B' aup = í/2(B a B' + B'M - l/2(B a B p - B' a B' ). (5) 

First, we prové that the following inequality proposed by 
Uffink for qubit systcms arc also vàlid for an arbitrary 
dimcnsional system: 



(B a ) 2 + (B' a ) 2 <2l Q l- 1 ,(|a| >2). 



(6) 



In order to see this, we use the following lemma. 

Lemma: Let — 1 < X{,X[ < 1 be Hermitian operators 
fi = 1,2 ), and define Y, Y' as follows, 



f(Y,Y') = f(X 1 ,X' 1 )®f(X 2 ,X' 2 ). 



(7) 



Then 



(Y) 2 + (Y 1 ) 2 < 2. (8) 

The lemma is proven in the following way: Note that 

Y = (1/2){Xí(X 2 + X' 2 ) + X[(X 2 - X' 2 )}, 
Y' = {1/2){X' 1 {X 2 + X 2 ) + X Í {X' 2 -X 2 )}, (9) 

and let Bg be Y cos9 + Y' s'm9. Let us derive the maxi- 
mum value of tr[pB§]. Note that tr[pBg] is a linear func- 
tion of each X^ or X[, keeping p fixed. Therefore we may 
consider only the set of extremal points in the convex set 
of Hermitian operators with — 1 < X < 1. Hence we 
may assume X 2 = X' 2 = l(i = 1,2). We thus have 



Y 2 = Y' 2 = 1 - (1/4) [X x , X{] ® [X 2 ,X 2 ] 
= l + Aï®A 2 , 



where A i = [i/2)[Xi,X' i ] (A i are Hermitian operators) 
and 

{Y, Y'} = (1/2){X X , X[} ® {X 2 , X' 2 } = 2A+ ® A+, (11) 

where = (1/2){Xí, X^} (Af are Hermitian opera- 
tors). Then we obtain 

Bg = 1 + Al <g> A 2 + sin 29 A^ ® A 2 . (12) 

This implies 

tr[/oSg] < max{l + \\Af <E> A^ + A^ <g> A 2 \\ , 

1 + \\A+ ®A+ - A^ <z>A 2 \\}, (13) 

where |j • || means the operator norm. Note that 

Af (g> A^ ± Al & A 2 
= (1/2){(A+ +iAl)®{A+TiA 2 ) 

+(A+ -iAl)®{A+±iA 2 )}, (14) 

and 

A^ (g) A\ + A{ <g> A 2 
= {l/2){X[Xi <g> X 2 X 2 + X X XÍ ® X£X 2 ). (15) 

According to relationships such as = 
\\(X 1 X' 1 )(X' 1 X 1 )\\ í / 2 = \\l\\ í / 2 = l,weget 

\\A+ ®A+ + A^ ®A 2 \\<1. (16) 

Similarly, we also get 

\\A~l ®A%- Al ® A 2 \\ 
= (1/2)\\X' 1 X 1 ®X 2 X 2 + X 1 X' 1 ®X 2 X 2 \\ < 1. (17) 

Therefore we have |tr[p_Be]| 2 < tr[pi?g] < 2 by the vari- 
ance inequality. Now by taking 



cosí 



(Y) 



i, Slïlf 



(Y>) 



y/(Y) 2 + (Y')< 



(18) 



we obtain (Y) 2 + {Y') 2 < 2, Q.E.D. 

Let us consider a set a C N„, where \a\ > 2. Let 7 be 
a\{j}, where j G a. Then, from Eq. (|j), we have 

f{B a ,B' a )^f{B 1 ,B' 1 )®f(A J ,A' J ). (19) 

It has been known that the maximum of (£> 7 ) is 2^' 1 ^ 1 ^ 2 
]Ï0f. Noting that f(cx,cy) — cf(x,y),c g R and —1 < 
S^ItI- 1 )/ 2 » < i j according to the lemma by taking 
Xl = 2-dTl- 1 )/ 2 ^,^ = 2~^\^I 2 B\,X 2 = A h and 
X 2 = A[j , we obtain the quadratic inequality 



(B a ) 2 + {B'J 2 < 2H- 1 , (H>2), 



(20) 



(10) 



where we used |a| = I7I + 1. 

Next, we calculate an upper bound of (Kn„) 2 + (B^ ) 2 
for states of the form l S)f— 1 p ai . From Eq. (^), we have 
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(B aU p) 2 + (B' aU/3 ) 2 

= (l/2)((B a B' + BM 

Using Eq. @, we obtain 



(B a Bp-B' a B' (j ) 2 ). (21) 



(l/2)«B Ql # 



N„\q 



ai D N„\ai 



+(B ai B Nn \ ai - B' ai B' Nn \ ai ) ) 
= (1/2){((6 Q1 )(^ W ) + (^X^w» 2 

+ ((6 Q1 )(S N , AQ1 )-(^ 1 )(^ W )) 2 } 
= (1/2)((6 Q1 > 2 + (^^((^w) 2 + ^N„\ Q1 > 2 )- (22) 
Repeating this, we obtain 

<#N„) 2 + <^nJ 2 = (V2)*- 1 ri«S Ql ) 2 + (^ ! ) 2 )· ( 23 ) 



Without loss of generality, we assume that \cti\ = 1 for 
1 < i < m and |ai| > 2 for m + 1 < i < k. Applying 
(B a y + (B' a ) 2 < 2 for \oti\ = 1 and Eq. (§), we obtain 



<Sn„) 2 + (^„) 2 

k 

< Yl 2 (l a 'l- 1 )(l/2) fe " m " 1 

i—m-\-l 



■\n-\-m— 2/c+l 



(24) 



where we used y^ =m+1 (|ai| — 1) = (n — m) — (k — m). We 
then conclude |Ï4| that, for any state p that is fc-separable 
with respect to a\, . . . , ctfc, 



(tr[pfí N J) 2 + (tr[pSi 



N„ 



< 2 



n+m-2fc+l 



(25) 



The maximum depends only on two parameters k and 
m but not on the detailed configuration of the partition. 
Cleary the bound ( ^H| ) is optimal. 

The inequality for testing full n-partite entanglement 
for n > 3 is obtained by maximizing the right-hand side 
of ( p5[ ) under the condition k > 2. Noting that m < k—1 
when fc < n, we obtain 



N„ 



> 2 < 2 n " 2 



(26) 

6|) imply full n-partite en- 



Violations of the inequality 
tanglement. 

For multiqubit systems, Ufhnk considered the case 
of partitions of the form {!}, {2, . . ■ , n\, and has pre- 
sented the quadratic inequality ( |26| ) for testing whether 
n-particle states are fully entangled ||. In what we 
should pay attention to, we have to check that for all 
hybrid separable-inseparable states except genuinc fully 
entangled states, the optimal upper bounds are smaller 
than or equal to 2™~ 2 , in order to ensure that the relation 
( p6| ) can be used as tests for full n-partite entanglement. 
In this point, we genuinely proved that the violations of 
the relation ( |2"6| ) are a sufficient for confirming fully n- 
partite entangled states. We have also proven that the 



relation (26) can be derived not only for multiqubit sys- 
tems but also for higher dimcnsional systems. 
The inequality (Eòl) also implies 



(27) 



It is known that |tr[jO,BN,J| < 1 when the system is fully 
separable p0[. Hence we obtain an upper bound p5| 



MpB 



N„ 



< 



) (n+m-2fe+l)/2 



1 



k < n 
k — n. 



(28) 



According to Eq. (||), the equality of the relation ( |2q ) 
holds when (B ai ) = (B' a .) = 1 for 1 < i < m, (B ai ) = 
(B' ai ) = 2d Q -l- 2 )/ 2 for m + 1 < i < k- 1, and (B ak ) = 
2(l a <=l- 1 )/ 2 . We can find a state and Hermitian operators 
— 1 < Aj,A'j < 1 that satisfy the above relations [^6) . 
Hence the bound (|2|) is optimal. 

For partitions of the form {!}, {2}, . . . , |m|, jm + 
1, . . . , n}(m < n — 1), the relation ( |2q ) leads to the result 
of Gisin and Bechmann-Pasquinucci i.e., the bound 
\{B-H n )\ < 2("- m " 1 )/ 2 . Noting that m < k - 1 when 
k < n, the relation ( |28| ) leads to the result of Werner 
and Wolf JÏÏj], i.e., |(B N „)| < 2("- fc )/ 2 by taking the 
maximum over m with fixed k. Collins et al. consid- 
ered the cases for partitions of the form {1}, {2}, {3, 4} 
or {1,2}, {3, 4} or {1}, {2, 3, 4} and presented the bounds 
as -\/2, V^, 2, respectivelyj 7J . These bounds are also de- 
rived from the relation (|28). 

So far, we derived the threshold value (i.e., 2™ -2 ) of 
(^n„) 2 + (B^ ) 2 for use as a full n-partite entanglement 
witness over all observables satisfying —1 < Aj,A'j < 1. 
Now, let us consider an additional assumption that the 
two measured observables anti-commute, i.e., {Aj, A'j} = 
OVj. This assumption is approximately fulfillcd within 
the accuracy of the measurement apparatus in the com- 
mon experimental situations, e.g., when we measure 
Pauli operators a x and a y for each particle. With this 
assumption, the threshold value of (2?n„) 2 + (^ísr,,) 2 be- 
comes as small as 2™~ 3 as is shown below. This implies 
that we can use a stronger quadratic inequality as tests 
for full n-partite entanglement in this case. 



Suppose that {Aj,A'A = and 



1 < A j, A j 



< l.Vj. 



Let us take Aj g — Aj cos 6 + A'jSxa.6, and derive the 
maximum of the vàlues tr[pA,- e ]. Since we are only inter- 

ested in the maximum, we may assume Aj 2 = A'j = 1. 

Then we get A/ g = 1 + (l/2){A. h A'j} sin2(9 

s to \tr\pA 

1. Now take cos 9 = (Aj) / ^{Aj 



1. The 

2l 



variance inequality leads to |tr[pA,- e ] | 2 < tr[pv4 Je ] = 



(4: 



sm f = 



(A'^/^A,) 2 + (A^)*, then we get (A,) 2 + (A^) 2 < lVj. 



This means that the relation (^) holds even for \a\ 
Hence we obtain 

<S Q ) 2 + <^) 2 <2H-\(H>1). 



(29) 
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Similar to the argument that derives 
relation (|29l) , we conclude 



50, applying the [1] 



(tr[pfí : 



(tr[p^J) 2 < 2 



n-2fe+l 



(30) 



The inequality for testing full n-partite entanglemcnt 
is obtained by maximizing the right-hand side of ( J30| ) 
under the condition k > 2. We obtain |ïjj 



N, 



N, 



< 2 



n — 3 



(31) 



We give an example that the relation (|T]) is stronger 
than (^) as a witness of full n-partite entanglement for 
multiqubit systems. We assume that A,- — dj ■ a, A' 



■ (T, whcrc d'i and 



are normalized vectors in R 3 and 



a is the vector of Pauli matrices, i.e., a = {<J x ,o'y,o'z)· 
The condition {Aj,A'j} = leads to dj -aL = 0. Let us 
consider the following multiqubit states JTÏJ : 



p = x\<ï> n )(<ï> r 



1 



2" 



(32) 



where / is the identity operator for the 2"-dimensional 
space and |<& n ) is an n-qubit GHZ state |Q, i.e., 



1$, 



1* 



+l,+2,...,+n) + |- 



,))■ (33) 



It is easy to see that the maximum of (£>n„) 2 + (£>n 



is 2™ 1 x 2 with d'i ■ a'- 

3 

that x is in the range of 



2 + (V 

OVj (See |19|). Hence, assuming 



1 1 

2 <X -V=2> 



(34) 



we can confirm full n-partite entanglement from (31), 
which cannot be confirmed by (|2(j). Hence if the mea- 
surement setups are precisely chosen as {Aj, A^} = OVj, 
then one can use a stronger inequality as tests for full 
n-partite entanglement in comparison with the relation 

In real experimental situations, we cannot claim that 
{A,,-, A^} = with arbitrary precision. The relevance of 
the bounds claiming full n-partite entanglement assum- 
ing that |({A.,-, A^-})| < e, where e means experimental 
errors, would be worth further investigations. 

In summary, we have derived the quadratic inequality 
that is utilizable to test full n-partite entanglement not 
only for qubit systems but also for higher dimensional 
systems. This helps the analysis of experimental data in 
realistic situations. We have also shown that when the 
two measured observables are assumed to precisely anti- 
commute, we can use a stronger quadratic inequality as 
a witness of full n-partite entanglemcnt in correlation 
experiments. 

K.N. thanks Masanao Ozawa for helpful discussions. 
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